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INTRODUCTION 


The  present  paper,  prepared  under  Navy  Contract  N6ori-201, 
Task  Order  No.  1,  is  the  third  of  the  series  announced  in  the 
report  Waves  Over  Beaches  of  Small  Slope.  Under  a  Dock.  Under 
an  Overhanging  Cliff,  and  Past  Plane  Barriers.  (IMll-NIU  l6l). 
Beyond  the  solution  of  the  specific  problem  indicated  in  the 
title  it  contains  new  observations  and  methods  which  may  be 
interesting  from  the  general  point  of  view  of  mathematical 
analysis. 

R.  Courant 


WAVES  IN  THE  PRESENCE  OF  AN  INCLINED  BARRIER 

By  Fritz  John 

The  present  paper  deals  with  the  motion  of  water 
of  infinite  depth  in  the  presence  of  a  fixed  barrier 
of  length  A  which  extends  downwards  from  the  siirface 
and  forms  an  angle  -kt   with  the  horizontal  direction 
where  n  is  a  positive  integer.  For  the  case  of  a 


vertical  barrier  (n  =  1)  the  motion  was  determined 
by  P.  Ursell  [1].  That  author  obtains  his  solution 
by  superposition  of  exponential  solutions,  arriving 
at  an  integral  equation  for  the  unknown  function  on 
the  vertical  line  containing  the  barrier.  He  succeeds 
in  solving  that  integral  equation  explicitly  and  ob- 
tains a  solution  which  is  characterized  by  the  fact 
that  the  velocity  of  the  liquid  becomes  infinite  only 
at  the  edge  of  the  barrier. 

The  methods  developed  by  J.  J,  Stoker  and  H,  Lewy 
for  determining  waves  on  sloping  beaches  (see 
IMF-NYtJ  161*)  will  be  seen  to  furnish  a  simpler,  and 
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in  some  respects  more  complete  solution  even  in 
the  case  n  =  1  and  can  be  extended  readily  to 
general  n.  By  the  use  of  complex  functions  the 
determination  of  the  motion  is  reduced  to  the  solu- 
tion of  an  ordinary  differential  equation  with  con- 
stant coefficients  instead  of  an  integral  equation. 
The  rSle  of  the  various  physical  assvimptions  in 
uniquely  determining  the  solution  is  brought  out 
clearly,  and  assumptions  that  have  to  be  made  in 
order  to  make  Fourier  integral  transformations  valid, 
are  avoided. 

The  method  pei*mits  us  to  find  the  most  general 
solution  with  singularities  of  a  prescribed  nature. 
Thus  we  shall  consider  solutions  that  in  addition  to 
having  infinite  velocity  at  the  edge  of  the  barrier, 
have  singularities  of  the  source  or  sink  type  at  the 
points  v/here  the  free  svirface  reaches  the  barrier. 
Following  the  example  of  J,  J.  Stoker  [2]  the  occ^ir- 
rence  of  such  singularities  may  be  thought  to  corres- 
pond to  the  breaking  of  waves  at  the  barrier  with 
corresponding  loss  of  energy  to  the  liquid.  Such 
breaking  is  likely  to  occxir,  where  the  barrier  forms 
a  small  angle  with  the  surface. 

The  motion  of  the  liquid  in  the  presence  of  a 
barrier  of  infinite  length  extending  downwards  from 
a  point  inside  the  liquid  can  be  treated  by  the  same 
general  method.  For  n  =  1  the  problem  has  already 
been  solved  by  Ursell  [1],  The  complex  fiinction 
method  permits  in  fact  to  solve  a  great  variety  of 
problems  for  water  waves  explicitly.   In  a  later 
paper  the  author  will  take  up  an  application  to  the 
problem  of  determining  the  motion  of  floating  bodies. 

The  solution  of  the  problem  concerned  with  plane 
barriers  in  the  general  case  will  be  carried  out  here 
to  the  point  where  the  problem  is  reduced  to  the 
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solutlon  of  a  system  of  5n-l  ordinary  linear  equa- 
tions for  parameters  c,  on  which  the  motion  depends. 
The  coefficients  of  these  equations  are  certain 
transcendental  functions  of  the  length  A  of  the 
barrier.   It  did  not  seem  feasible  to  carry  out  the 
solution  of  these  linear  equations  explicitly  for 
general  n. 

In  the  case  of  a  vertical  barrier  (n  =  1)  the 
solution  of  these  linear  equations,  whose  coefficients 
are  now  Bessel  functions,  is  quite  simple  and  is 
given  here.   The  case  of  a  barrier,  forming  a  45° 
angle  with  the  surface  (n  =  2)  has  been  carried  out 
completely  by  E.  Isaacson  and  B.  Grossman,  The 
resiilts  of  their  reduction  are  given  in  a  section 
elsewhere.  The  coefficients  here  are  again  inte- 
grals of  products  of  exponential  and  algebraic 
functions,  and  may  be  considered  as  generalizations 
of  Bessel  functions. 

The  author  is  greatly  indebted  to  Professor 
J.  J.  Stoker,  whose  work  forms  the  basis  of  the 
present  investigation,  and  whose  advice  during  its 
progress,  particularly  relating  to  the  physical 
interpretation  of  the  results,  has  been  indispensable. 
The  results  and  notations  of  the  preceding  report 
IMM-NYU  161  will  be  constantly  referred  to  in  the 
course  of  this  paper. 
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1,  Method  for  finding  the  general  solution. 

We  adopt  the  notation  of  the  IMM-NYU  report 
n\amber  161  and  use  the  diniensionless  variables 
introduced  there  in  (2,17),  which  again  will  be 
denoted  simply  by  x,y,t.  In  those  variables  the 
barrier  extends  from  the  origin  to  the  point 
(x-j_,yj_),  where 

(1.1)  ^1  ~  ^  °°^  ^  »  y  =  -a  sin  ^  . 

Here 

(1.2)  a  =  mA  =  -^  A 

where  A  is  the  length  of  the  barrier  and  cr  the  circu- 
lar frequency  of  the  motion  in  the  original  variables. 

The  velocity  potential  $  (x,y,t)  is  a  poten- 
tial fxinction  defined  in  the  region  T*  bounded  by 
the  X-axis  and  the  barrier.  Besides  satisfying  the 
potential  equation 

(1.3)  $xx  ^  lyy  =  0  , 

$  is  assumed  to  be  simple  harmonic  in  the  time: 

(1.4)  i^^  +    $  =  0  . 

It  is  assxuned  that  $  and  its  first  partial  deriva- 
tives are  continuous  on  the  boxjndary  of  T'  ,  with  the 
possible  exception  of  the  edge  (x,,y, )  of  the  barrier 
and  at  the  points  (+  0,0)  where  the  free  surface 
reaches  the  barrier.   (In  what  follows  these  three 
points  will  be  referred  to  as  "comers"  of  T^  ,) 
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On  the  free  surface  7=0  the  condition 


(1.5) 


-  I  +   §y  =  0 


is  imposed.  Along  the  barrier  the  normal  velocity 
shall  vanish: 


(1.6) 


^n  =  0 


The  motion  still  depends  on  the  assxmiptions  mad© 
on  the  behavior  at  infinity  and  at  the  three  exception- 
al corner  points,  llies©  assiimptions  will  be  expressed 
here  as  conditions  on  the  velocity,  which  has  direct 
physical  meaning,  instead  of  conditions  for  $  itself. 

It  will  be  assumed  that  the  components   ^ 
and   $  stay  bounded  at  infinity  and  tend  towards  0, 
as  y  ->  00  , 

At  the  two  comers  of  "T  at  the  origin 


(1.7) 


(x^  +  y^)($ 


X   *y' 


shall  stay  botmded,  and  at  the  edge  (x,,y,) 
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:i.8)         lim      (U-x^)^  +  (7-7i)^(  1^  +  fy)  =  0  . 
(x,y)->(Xj_,y3_) 

Condition  (1,4)  implies  that  $  is  of  the  form 

(l.Sa)   $(x,y,t)  =  <P-L(x,y)cos  t  +  <p2(x,y)sin  t   . 

Here  the   <Pv(x>7)  are  potential  fxmctions  of  x,y  in 
satisfying  the  same  regularity  conditions  and  bo-undary 
conditions  as  $  .  Formula  (l.Sa)  corresponds  to  a 
representation  of  ^  as  the  sum  of  two  "standing 
waves".   It  is  obvious  that  the  potential  functions 
(p,  and   (Pg  separately  satisfy  the  bOTxndary  condi- 
tions (1,5),  (1,6)  and  the  assumptions  on  the  behavior 
of  $  at  infinity  and  at  the  singular  points. 

It  will  be  shown  that  a  potential  fimction  <p 
defined  in  "P  ,  that  satisfies  the  free  surface  con- 
ditions (1,5)  and  for  which  (p      and  (P     are  bounded 

X       *  y 

at  infinity  and  vanish  for  y  -*•  -co  ,  is  of  the  form 

(p(x,y)  =  e^(oC  cos  x  +  p  sin  x)  +  R 

where  R  and  its  derivatives  tend  to  0  for  (x,y)  ->  od  . 
Here   «<•  ,3  are  constant  which  may  however  be  differ- 
ent for  positive  and  negative  x.   Correspondingly 
^ (x,y, t)  for  large  positive  x  is  asymptotically 
represented  by  an  expression 

e^[A  cos(x+t)  +  B  sin(x+t)  +  C  co3(x-t)  +  D  sin(x-t)] 

and  for  large  negative  x  by  an  expression 

e^[A«cos(x+t)  +  B'sin(x+t)  +  C»cos(x-t)  +  D«3in(x-t)]  . 
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Thus  ^  can  be  considered  approximately  for  large 
Ix|  as  a  svan.   of  two  progressing  waves* 

The  values  of  A,B,C,D,A',BSCSI>'  are  not  in- 
dependent from  each  other  but  are  linked  together 
with  the  singularities  of  $  through  the  energy 
integral.  We  have  indeed  the  following  theorem: 

ir  <|  satisfies  conditions  (1.3),  (1.5),  (1.6) 
and  if 


II 'f 


it 


)dxdy 


exists  over  any  finite  subregion  of  I  ,  then  the 
values  of  A,B,C',D*»  already  uniquely  determine 
C,D,A',D'.   That  is,  the  incoming  waves  uniquely  de- 
termine the  outgoing  waves  at  infinity  if  no  sources 
or  sinks  of  energy  exist  in  P  . 

For  the  proof  of  this  statement  it  is  siifficient 
to  consider  the  case  where  A,B,C',D'  vanish.  Con- 
sider the  slit  rectangxolar  region  y,  whose  boundary 
b  consists  of  the  vertical  line  segments  x  =  x^, 
X  =  Xr,,  the  horizontal  line  segment  y  =  y  ,  the  free 


X-SCg 


y^Vo 


surface  and  the  barrier.  The  expression 


-8- 


2 


y=0 


suggests  itself  as  the  total  amoiont  of  energy  con- 
tained in  Y»  ill  view  of  the  fact  that   $^  repre- 
sents the  surface  elevation  in  the  units  chosen  here. 
We  have 

i  =  //($x  §xt  -^  #y  ^yt)^'^^  ^    j     ^  ^tht^ 

y=o 


(  ^  H  ^-  -  / '  *t  ^tt  ^ 


^1 

y=0 


-^2 


=' /  ^t  $x  ^y  -  /  ^t  ^x  <^y  -^  f  ^  ^y  -^  *tt)^  • 

/  7o  ^o  ^1 

x=Xp  x=x,  y=0 

Here  the  last  integral  vanishes  because  of  (1,4); 
for  large  positive  x-,  and  large  negative  Xg  the  first 
two  integrals  becorae 

-(1-e  °)[(Csin(x2-t)-Dcos(x2-tJ)^+(A«sin(x^+t)-B»cos(x^+t))^]  +  e 

where  the  error  £  involved  tends  to  0  as  x-, ->-od  , 
Xp-^'+CD  •  Then,  as  I  is  obviously  a  periodic  function 
of  t,  as  $  Is  periodic. 
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2Tr 
0  =  I(2Tr)  -  1(0)  =   f   ^  dt 


2lT 

Tr[A'^  +  B«^  +  C^  +  D^]  +   r   edt 

o 


As  £.  Is  arbitrarily  small  for  a  suitable  choice  of 
^1**^2'  ^*  follows  that  A*,B',C,D  vanish. 

At  the  end  of  this  paper  it  will  be  proved  -chat  a 
solution  (p  which  tends  to  0  for  x  -^  oo  and  x  — >  -co  , 
has  to  vanish  everywhere.   It  then  follows  that  the 
motion  of  the  water  in  the  presence  of  the  barrier 
is  uniquely  determined  by  waves  coming  in  from  infinity 
on  either  side,  as  long  as  no  "logarithmic"  or  higher 
singularities  are  present.  More  generally  the  motion 
will  be  determined  by  the  incoming  waves  and  the 
"strength"  of  the  energy  sources  and  sinlcs  at  places 

-v. 

where  they  exist. 

In  accordance  with  the  method  outlined  in  the 
Bffl-NYU  T'eport  number  161  the  above  conditions  for  the 
<p,  (x,y)  are  translated  into  conditions  for  the 
analytic  function  "X  (z)  =  "X.  (x  +  iy),  whose  real 
part  is   9t^,  X  is  determined  uniquely  by  (P^ 
except  for  an  arbitrary  pure  imaginary  constant. 
Then  X  (z)  is  regular  analytic  in  the  interior 
points  of  F  ,  "X'Cz)*  and  hence  also  %(z),   is 
continuous  on  the  boundary  of  T  with  the  possible 
exception  of  the  origin  and  the  point  at  =  x-  +  iy-j_, 
where  we  put 


"This  proof  is  comparatively  involved  and  makes  use 
of  the  special  nature  of  the  region  P  •  It  is 
plausible  that  an  analagous  theorem  exists  for  the 
motion  of  a  liquid  in  the  presence  of  any  type  of 
obstacle. 
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•rri 
(1.9)  £=e-^ 


The  assiimptions  made  on  the  velocity  guarantee 
that 

(Al)       X'(2^    ^^  bo\mded  as   z    ->  oo    in    F  and  that 

lim     X*^^  +  iy)   =  0 
y  -5»-Q0 

(A2)   that  z  X'(2i)  is  boimded,  as  z  ->  0  in  V   and 
that 

(A3)    (z  -  afc)  "X' (z)  ->  C  as  z  ->  at  in  r 

The  boundary  conditions  become,  in  accordance  with 
the  report  mentioned  above, 

(1.10)  Im(^  +  i)  X(z)  =0  ^o^  ^^^1  ^  ^  ^ 

(1.11)  IviitXHz))   =  C  along  both  sides  of  the 

barrier,  i.e.  for  z  =  tfc  , 
where  t  is  real  between 
0  and  a. 

The  essential  step  in  finding  a  fxinction  %,  {z) 
satisfying  all  the  conditions  mentioned  consists 
in  reducing  the  two  boundary  conditions  (1.10), 

(1.11)  to  identical  form  by  introducing  the  auxiliary 
function 

(1.12)  P(z)  =  D''(D-C3^)(D-f2)...(D-fn)X(z)  =  L(D)  X(z) 
Here  D  denotes  the  differential  operator 
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D  = 


^ 


and  the  constants   ^  v  ^.re  defined  by 

(1.13)        Ck  =  ^^    ~  ^® 

We  notice  the  following  properties  of  the   4*  ij.> 
which  are  essential  for  all  that  follows: 


(1.14) 


R(4ij.)  ^  0  i'or  k  =  0,1, ...,n 


(1.15) 


^n  =  -i   '   ^o  =  ^ 


(1.16) 


^k  "  ^n-k 


(1.17) 


^k=  ^^n-: 


k+l 


The  symmetry  conditions  (1.16),  (1.17)  imply  iiie 
following  properties  for  the  polynomial  L(D): 


n-1 


(1.18) 


P(D)  =  ^  =  D^  It  (D-^i^) 
""'^       k=l     ^ 


D-f 


n 
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Is  a  polynomial  with  real  coefficients,  as  the  roots 
are  conjugate  imaginary  in  pairs.  Similarly 

n 

(1.19)  L(D)  =  -(£D)^  TT  (£D-e?j^)  =  Q(£D)  •  (6D) 

k=l 

where  Q  is  a  polynomial  with  real  coefficients. 

As   7(  ♦  +  i  X  is  real  and  continuous  along 
the  real  axis,  it  is  regular  there  as  well.  Then 

(1.20)  F(z)  =  P(D)  (  X»  +  iX) 

is  regular  along  the  real  axis  and  also  real,  as 
P(li)  is  a  polynomial  with  real  coefficients.  Similar' 
■^7  ^X'  ^2  real  and  continuous  along  each  side  of 
the  barrier,  and  hence  regular  there,  though  the 
values  on  one  side  will  not  be  the  analytic  continu- 
ation of  those  on  the  other  side.   It  follows  again 
that 

(1.21)  P(z)  =Q(6D)(e.X')  =  Q,i-^)itxHtt)) 

is  regular  and  real  along  the  barrier. 

Thus  the  function  P(z)  is  real  along  the  whole 
boxandary  of  T  ,  with  possible  singularities  at 
0,00, a£.   It  remains  to  discuss  the  behavior  of  P 
at  those  points.  The  function 

(1.22)  f(z)  =  (D  +  i)  X(z) 

being  real  along  the  real  axis,  can  be  continued 
across  the  axis  by  the  forr:ula 

(1.23)  f (z)=  fTzl 
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and  becomes  regular  in  the  region  I  '"  formed  by 
adding  to  V   the  points  obtained  by  reflection  with 


respect  to  the  real  axis.  By  assumption  (Al)  there 
exists  a  constant  M  such  that 


for  large  z  in  "F  ,  Then 


|f(z)|  =  I  X»  +  Itl  <  M|z|  +  C 

for  large  z  In  F"  ,  and  hence,  because  of  (1,23), 
in  a  fiill  neighborhood  of  infinity,"  Thus  f (z)/z 
is  bounded  at  infinity,  and  hence  regular  there: 

-2 

f(z)  =  d  ,z  +  d^  +  d,z  +  dgZ   +  •••   • 

As  X*  ^^  bounded  at  infinity  in  T^  ,  it  follov/s 
that  in  V 


""'This  conclusion  would  have  to  be  modified  for  the 
case  of  a  submerged  barrier  extending  to  infinity. 


-14- 


^(z)  =  -id  -iZ  +  (bounded  terms)  . 

Now  the  assumption  11m   7^'  (x  +  ly)  =0  implies 
that  ^  *-°° 

Hence  we  must  have 


d,^  =  0  . 


The  coefficients  d,  in  the  expansion  of  f  are  real, 

as  f  is  real  along  the  real  axis.  By  changing 

■J(^(z)  by  a  suitable  pure  Imaginary  constant  we  can 

bring  about  that  d^  vanishes*   Then 
°  o 

oo_      , 
(1.24)   f(z)  =  (D+i)X  =  ^  dj^z    for  |z|  >  a  . 


It  follows  that  F(z)  =  P(D)  f  is  of  the  form 

(1.25)      P(z)  =  >>    -VZ    for  hi  >  a 

k=n+l  *^ 

with  real  at,  ,  as  P(D)  is  a  real  polynomial  contain- 
ing the  factor  D  • 

To  study  the  behavior  of  F  at  the  comers  of 
T  ,  we  make  the  general  remark  that  if  g(z)  is 
defined  in  a  region  y  whose  boundary  contains  the 
origin,  and  if 

lg(z)|  £  jly        in  Y 
then  by  Cauchy»3  formula 
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Is""(-)l  g^ 


in  any  point  of  Y  that  has  a  distance  >  p  from  the 
boiondary  of  y» 

We  nov/  apply  this  argvirnent  to  the  right  hand 
corner  of  V   at  the  origin.  We  have  from  (A2) 


-^t  =  o(-)  near  0  in  T 


Hence  also 


X(z)  =  0(i)   and  f  =    ^    +  ±X  =   0(|)  . 


Then  also  from  (1.23) 


f  =  0(i)  near  0  in  T  '"" 


It  follows  that 


P(z)  =  P{D)f  =  O(-i-)   , 


2n' 


if  p  is  the  minimxmi  distance  of  z  from  the  boundary 
of  T* '",  i.e.  from  the  rays 


arg  z  =  -^  and  arg  z  =  - 


IT 
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Simllarly 

in  the  right  hand  corner  of  T  and  hence,  if   tX' 
is  continued  across  the  barrier,  also  in  the  sector 

"  "^  =  arg  z  ^  0 


Then 

P  =  Q(tD)  tV   =  o(-?H^ 

p 

for  the  points  of  that  sector  having  a  distance 

>  p  from  the  boundary.  It  follows  that 

P(z)  =  0(-^) 

for  all  points  of  T^  having  a  minimum  distance 

>  p/sin  -^  from  the  origin.   Thus 

z^^  P(z) 

is  bovmded  near  the  right  hand  corner  of  T^  ,  The 
same  argument  shows  that  P  has  the  same  behavior 
near  the  left  hand  corner.  Furthermore,  the  same 
argument  shovfs,  in  view  of  (A3),  that 

(1.26)         lim   (z  -  ae)   P(z)  =  0  . 
z  ->ae 

The  actual  determination  of   the  function  P 
will  be  achieved  by  mapping  the  region   P   of  the 
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z-plane  conformally  into  the  sector  T 


■rr 


0  <  arg  ^  <  -^ 


of   the  w-plane,   by  means   of  a  fvmction 


(1.27) 


z   =  K(w)    . 


LJ 


Here  we  can  still  prescribe  arbitrarily  the  images 
of  three  bo-undary  points.  We  postulate  that: 

w  =  0  corresponds  to  z  =  -0  (left  hand  corner  of  i  ) 
w  =00  corresponds  to  z  =  +0  (right  hand  corner  of  T"") 
w  =  t,"     corresponds  to  z  =  a£  (edge  of  barrier). 

The  free  siu?face  in  the  z-plane  then  goes  over  into 
the  positive  w-axis,  the  barrier  into  the  ray 
arg  w  =  -gr  •  The  point  z  =  oo  goes  over  into  some 
point  ytu  on  the  positive  real  axis.  Prom  the  ch-ange 
of  angles  at  the  vertices  of  T""  we  conclude  that 
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yrr    \  2n-l 

near  w  =  0  :   K.(w)  ^  w 

near  w  =  oo  :   K(w)  ^  w" 

near  w  =6   :  K(w)-at  -  (w-fc  ) 

near  w  =/«'  :   K(w)  -^  (w-//)~ 

By  the  mapping  (1,27)  P(z)  goes  over  into  a 
fvinctlon 

(1.28)  G(w)  =  P(K(w)) 


v/hich  Is  regular  in  T   ,   and  real  on  the  boundary  of 

-1 

T  ,   with  possible  singularities  at  w  =  0,oo,  £   and 

ij,»     The  region  T  now  has  the  special  property  that 
a  function  G,  which  is  regular  in  T^  and  real  on  its 
boundary,  is  univalued  in  the  whole  w-plane.  G(v/) 
can  be  continued  into  the  whole  w-plane.  Reflection 
on  the  lines  arg  w  =  -^  and  arg  w  =  0  shows  that 

G(e^w)  =  G(w)=  G(w)  . 


2n 
Consequently  G  is  a  univalued  function  of  w   •   The 

only  possible  singularities  of  G(w)  in  the  whole 

plane  are  the  points 

w  =  0 

w  =  00 

^  -  A   -  c -l-2r 
w  =  A^  =  e 

-2t> 
w  =  B^  =  £  "> 
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whlch  must  be  isolated  singularities  and  not 
branchpoints.  As  z   F(z)  was  seen  to  be  bovuided 
near  the  origin  of  the  z-plane,  the  function  G(w) 
can  only  have  a  pole  of  order  2n(2n-l)  at  w  =  0 
and  a  pole  of  order  2n  at  w  =  oo  •  Because  of  (1,26) 
G(w)  can  only  have  a  pole  of  order  <  4n  at  w  =  6 
and  hence  at  any  point  A^,  Because  of  (1.25)  P(z) 
has  a  zero  of  order  n  +  latz=ao,  and  hence  G(w) 
has  a  zero  of  order  n  +  latw=/i.,  and  hence  at 
all  B  .  As  G(w)  has  poles  as  only  singularities,  it 

must  be  a  rational  function.  As  G  depends  only  on 

2n 

w  ,  it  must  be  of  the  form 

^+1         5n^    ,2nk 


C  2n  2n-,n+l        5n- 


V 


where  the  coefficients  Cj^  are  real,  as  G  is  real 
for  positive  real  v. 

The  mapping  fionction  K(w)  is  easily  determined. 
The  expression  z.w  =  wK(w)  is  obviously  real  on  the 

boundary  of  T_.  Thus  wK(w)  is  a  univalued  function 

2n 
of  w  for  all  w,  depending  on  w   alone.   In  addition 

K  has  a  zero  of  order  2n-l  for  w  =  0,  a  pole  of  order 

1  at  zx.  and  hence  at  every  B,^,  and  has  a  simple  zero 

at  w  =  00 .  Consequently  K  is  of  the  form 

with  a  constant  C.  This  constant  is  detei-mined  by 
the  condition  that  K(e  )  =  afc  : 

C  =  a(l  +/Ji?'^)    . 
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Pinally  the  constant  yw.  is  now  seen  from  the 
condition 


K»(e"-)  =  0 


to  have  the  value 

1 
(1.30)  />o=  (2n-l)"  "^ 


Thus  the  mapping  is  given  by 

2n-l 

(1.31)         z  =  K(w)  =  ^-^^^ -g.--   . 

(2n-l)w^^-l 


Thus  P  is  of  the  form  (1,29)  with  jju  given  by 
(1,30),  and  w  and  z  being  connected  by  the  algebraic 
relation  (1,31),  F  is  seen  to  depend  on  5n-2 
constants  c,^.  The  exact  determination  of  P  and  hence 
of  X  then  only  requires  a  closer  analysis  of  the 
resulting  expression^  and  the  derivation  of  the 
conditions  imposed  on  the  c,,  by  the  original 
assumptions  on  X  •  We  will  first  carry  out  this 
analysis  in  detail  for  the  case  of  a  vertical  barrier 
(n  =  1),  and  less  completely  for  the  general  case 
in  the  next  sections. 
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2,  The  case  of  a  vertical  barrier, 


In  the  case  n  =  1  formula  (1,12)  reduces  to 


(2.1) 


P(z)  =  D(D  +  1)  7[(z)  . 


con- 
sisting of  the  z-plane  slit  from  -ia  to  +ia,  P  is 
real  along  the  real  axis  and  along  both  sides  of  the 
nlit.  The  mapping  function  K  becomes  here  according 
to  (1.31} 

(2.2)  z  =%^  =  K(w)  . 

w  -1 


This  function  maps  the  first  quadrant  T  of  the 
w-plane  on  T^  ,  and  the  right  half  of  the  w-plane 


on 


".  Solving  for  w  we  get  for  z  in  I 


(2.3) 


w  =  ■  ■»■. 


a^+z^ 


/~S — S'  T-«  -J'- 

where  V^-  +2   denotes  that  branch  regular  in  l  " 

that  behaves  like  z  at  infinity. 
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Tiie  fiHiction  G(w)  becomes  here 
(2.4)      G(w)  =  V-       -l  -Tf    (c^+  c,w  +  CpW  +  c,w  )  . 

Now  from  (1.11)  1  "Viz)    is  real  along  the  barrier, 
and  hence,  as  a  function  of  w,  is  real  along  the 
positive  iraaginary  w-axiS'.  As  by  assxamption  (A3) 

lim  (z  +  ia)  i  X'(z)  =0 
z-i»'-ia 

it  follows  that  i  X*  as  a  fimction  of  w  has  at  most 
a  pole  of  order  1  at  the  point  v;  =  1.   Then 

i  7('  =  -^^  +  (regular  function  of  w-i)  • 


Consequently 


-23- 


2 
nf„\   dz  _  -  2a (w  +1)  _  dz  ,  d  ^  ^n  ^, 


(2.4a)  =  ^  -x»  .  1  ^.  dz 


:     M  +  (regular  function  of  w-i) 
(w-i)'^ 


Thus  the  residue  of  the  function 


^^""^  ^  ^   "  w^(l4-W^)^  ^''°''  ''l''^"'  ''S^^''  °3''^^ 


at  the  point  w  =  i  must  vanish.   Consequently  we 
get  the  condition 


(2,5)         3c^j  -  c-j^  -  Cg  +  3Cg  =  0  . 


Substituting  the  expression  (2,3)  for  w  into 
(2.4)  we  obtain 

2                                          a(-3c  -c,+Cp+3c,) 
P(z)    =  V  [(c   +c,+C2+C3)    +  o     ^     2 3_ 

I/a  +z 
a2(  30^-0^^-02+303)        a^( -0^+03^-03+03) 

+     R „ +     J  , 

With  the  help  of  (2,5)  this  expression  may  be  written 

2c^a2(  lK?-a)        2c,a2(\P+?+a) 

F(z)   =       Q  +       ^  

2\/~S     S  2\r2     2 

z  Vz  +a  z   yz  +a 

(2.6)  3 

a  (-C0+C2-2C3) 


^ 


^ — S  3 
+a 
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Having  determined  F,  the  function   X(z)  is 
a  solution  of  the  differential  equation  (2,1). 
The  general  solution  of  that  equation  can  be  written 


1-1 


-1 


,-li 


X(z)  =  [D(D+i)]"-^  P  =  -i[D"-^  -  (D+i)  '  ]P 


=   -i(  I    F{^)dyi+^)    +  ie"^^(  I      e^^   P(^)'iV^  +  p) 


with  arbitrary  conplex  constants   o^^p^z^.  As  P(y|) 
behaves  like  Yl"^  at  infinity,  the  integrals  con- 
verge, if  we  choose  z  =  -oo  ,  and  prescribe  that  the 
path  of  integration  shall  lie  in  I       and  shall  approach 
-00  along  the  negative  real  axis.  Then 


/ 


(2,7) 


'X«(z)  =  e"'^[  f    e^''^  P(vj)dV|  +  p] 


"•CO 


(2.3)        -;{^   +  IX    =  P(y()di^  +  (K 

<^-oo 
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According  to  (1.24)  7[»  +  iX   must  vanish  at 
InTinity,   This  implies  obviously 

(2,9)  o6  =  0  . 

In  view  of  the  fact  that  F(y|)  =  ^^r^)*   ^^®  condition 

(2.9)  is  also  sufficient  to  guarantee  that 

5(_t  +  iX.  is  real  along  the  real  axis  and  vanishes 
at  infinity. 

Prom  the  lemma  that  will  be  proved  below  (see 
p.  45)  it  follows  that  the  expression   X*  given 
by  (2,7)  satisfies  in  J'' for  large  z 

X»  =  Oie'^^)    +  0(i)  . 

Hence  the  velocity  is  bounded  at  infinity  and 
vanishes  for  y  ->  -oo  •  The  only  condition  that 
remains  to  be  checked  is  the  reality  of  iX' 

at  the  barrier, 

— iz 
As  e    is  real  on  the  barrier,  we  must  have 

z 

(2.10)  R[  r  e^^P(>rj)dy|  +  p]  =  0 

''-00 


for  z  on  the  barrier  on  either  side. 

We  consider  a  point  z  =  it  with  -a  <  t  <  0 
on  the  left  hand  border  of  the  barrier.  The  lower 
limit  in  the  integral  in  (2.10)  may  be  replaced  by 
+ioo ,  as  long  as  the  path  of  integration  does  not 
cross  the  positive  real  axis.  Applying  integration 
by  parts  and  omitting  terms  with  vanishing  real  pai'ts 
we  find  from  (2,10),  (2.6) 
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R(P)    =  R(i 


it 


.^^[- 


4c_a 


+  21 


(0^-03)   a  1^ 


^+a^+a 


a(-c   +Cp-2c,)Yl 


\/ri 


^t;^^ 


2 


/ta 


ta 


Here   the  only  contribution  to 

it 


R(i   J       e^n    y|-l    dVI 


+ioo 


arises  from  the  origin,  which  contributes  the  value 
-TT.  The  remaining  terms  of  the  integrand  are  inte- 
grable  along  the  imaginary  axis.,  Taking  the  path  of 
integration  along  that  axis  we  notice  that  the  only 
contributions  arise  from  the  points  between  +ia  and 
ioo ,  and  only  from  the  terms  containing  a  radical. 
Substituting  Yl  =  i^  a  we  finally  obtain 
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.2  .  .,_   .  ^-2y^e-^^^l 


R(p)  =  A:vc^ar   +  2{c^'-c^)a.^   I     e  '^  tSi^  d^ 

'l 


(2.11)  00 


+  a2(-c^+Co-2c,)  /  e-^^  -^  d^T 

1        iP^ 


The  integrals  occurring  in  {2,11)  can  be 
expressed  in  terms  of  Bessel  functions  K  of  the 
second  kind,""  We  have  indeed 

/P         r"5 — 


a 


KQ(b)db 


00 

(2.15)        S2(a)    =    j      e-*^    —^  d^     =  IC^(a) 


•1  l/pl] 

In  this  notation 


(2.14)   R(p)   =  2Tra2(c^+C3)    +  2(cQ-C3)a^S3_(a) 


+   (-CQ+C2-2c3)a^S2(a)    • 


This  condition  then  is  equivalent  to  reality  of 
iX.*  on  the  left  hand  side  of  the  barrier.  The 
reality  of  iX*  on   the  right  hand  side  is  then  satis- 
fied as  well  as  a  consequence  of  relation  (2,5), 


""See  V/hittalcer  and  Watson:   "Modem  Analysis", 
p.  384. 
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Indeed  G(w)K»(w)  was  seen  to  be  of  the  form 
(2,4a)  near  w  =  i.  As  K»(i)  =  0,  the  expansion 
of  6^^^^'  in  terTns  of  (v-i)  contains  no  first 
order  term.  Then  the  function 


e 


^^("^^  G(v)K«(v) 


has  residue  0  at  v  =  i.  Then  the  expression 
z  w 

(3  +  r  e^"^  P(i\)dYl_  ""  P  "^  /  Q^^^^^  G(v)Kt(v)dv 

"^-00  1 


has  no  branchpoint  for  w  =  i.  Hence,  if  this 
expression  is  pure  imaginary  for  z  on  the  left  hand 
border  of  the  barrier,  i.e.  for  pure  imaginary  w 
between  0  and  +i,  then  it  is  pure  imaginary  for  all 
w  on  the  positive  imaginary  axis, 

(2.15)  lim  [X(z)  -  ipe"^'']  =  0  , 
z  -s»-oo 

To  investigate  the  behavior  of  X  (z)  for 
z  ->  +00  we  make  use  of  the  lemma  to  be  proved 
below  (see  p,  45  )•  We  obtain  that 

(2.16)  lim  [X(z)  -  i(!3  +  G)e"^^]  =  0 

Z-i>+00 

where 

(2.17)  C  =  dy  eM  F{^)    dY( 

the  integral  being  extended  over  a  closed  path  in- 
cluding the  slit  from  -ia  to  +ia.   It  is  permissible 
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to  obtain  the  asymptotic  behavior  of  the  derivatives 
of  X.  (2)  by  differentiation  of  the  formulae  (2,14), 
(2.13). 

We  have  then  the  result  that  the  most  general 
f "unction  %   (z)  satisfying  our  conditions  is  repre- 
sentable  in  the  form  (2.5a),  where  P(z)  is  given  by 
(2.5)  with  arbitrary  real  constants  Cq'^p'^S*  ^^^ 
where  z  =  -00  *  «<  =  0,  and  p  is  subject  to  the  con- 
dition (2.11),   The  behavior  of  this  solution  at 
infinity  is  then  determined  by  relations  (2.13), 
(2.14). 

Integrating  by  parts  in  (2.7)  we  obtain 


C  = 


-i^e^^[-2a2(cQ+C3)  i  +  2a(c^-C3) 


ff3 


+  a(-c^+C2-2c3)  —^  ]  dV^ 

=  -4Tra2(c^+C3)  -  2ai(c^-C3)  (^  (e^^-1)  ' ^  ^^  dY| 
-  a(-c^+C2-2c3)i(^  e^^  -pJl^  dY[ 


as 


/ 


M~+a 


dr|  =  0 


The  path  in  the  remaining  integrals  can  be  contracted 
into  the  slit  from  -ia  to  +ia,  and  we  obtain  for 
y|  =  ia^  ,  in  analogy  to  (2.14) 
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(2.18)   C  =  -4iTa^(cQ+C3)-4a%(cQ-C3)T-j^(a)-2ia^(-CQ+C2-2c3)T2(a) 


Here  we  put 

T 
(2.19) 


1         I "o  ^ 

l(a)  =  r  (e-^-1)  \ii^  d^  =  -TrlJ-L(ia)--rT  f   jQ(ib)db 
w-1  *'o 


1 


T^(a)  =  f   e"^  —i—  d^  =  -rrlJ,  (la)   . 

^       J-1      \fi:p 


The  general  solution  y^  (z)  is  seen  to  depend 
on  four  real  constants  c  jCgiC^  and  Iiii(p ) ,   It  is 
convenient  for  purposes  of  discussion  to  break  up  the 
solution  into  fo\ir  "fimdamental"  solutions,  each  with 
simple  characteristic  properties.  We  write 

Kz)  =  Im(p):;(Q(z)  +  4a2c3Xi(z)  +  4a%  X2(z) 

(2.20)  2 

+  2a  (CQ-C2+2C3)  A3(z)  . 

Here  each  fundamental  solution  "X  4 (z)  is  character- 
ized by  its  values  of  Im(p),c  .Cg^c^  as  follows: 

'Xq(z):   lm(p)  =1  ,   c^  =  Cg  =  C3  =  0 

"i  11 

A-l(z):   Im(p)  =0  ,  Cq  =  0,  Cg  =  — 5  ,  C3  =  — ^ 

^^(z):   Im(p)  =0  ,   Cq  =  -^,   C2  =  -ig,  ^3  =  0 

'X3(z):   Im(p)  =0  ,   c^  =  0,  Cg  = ig,  C3  =  0  . 
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Correspondlng  to  (2,20)  (2?  (x,y)  =R(7C(2))  can  be 
written  as  a  conbination  of  four  "fundamental" 
potential  functions  (pAx,j)t   which  are  the  real 
parts  of  the  liAz),     The  definitions  and  charac- 
teristic properties  of  these  functions  are  summarized 
on  p,  32 •  The  general  velocity  potential  as  a 
function  of  x,y,t  can  then  be  written  in  dimension- 
less  variables  as  a  linear  combination  of  the  (p-^ 
as  follows 

3 

(2.21)  $(x,y,t)  =]^  (fiC^cos   t  +  pj^  sin  t)  <j!>^(x,y) 

with  constants  «^ij.*Pv»  We  now  try  to  determine  the 
constants  in  such  a  way  that  the  resulting  motion 
corresponds  to  that  caused  by  a  progressing  wave 
arriving  from  x  =  -co  •  It  is  plausible  that  this 
wave  will  be  partly  reflected  by  the  barrier,  and 
partly  transmitted  to  x  =  +00  in  such  a  way  that  no 
"breaking"  of  waves  occurs  on  the  rigjht  hand  side  of 
the  barrier  and  no  wave  progressing  to  the  left  exists 
at  X  =  +00  • 

The  absence  of  breaking  waves  on  the  right  hand 
side  of  the  barrier  corresponds  to  the  assxjiirption 

^1  =  Pi  =  0 

as  the  term  with   <p,  is  the  only  one  ftirnishing  a 
logarithmic  singularity  there. 

For  large  negative  x  and  y  =  0  the  function  <jp 
is  represented  asymptotically  by  an  expression  of  the 
form 

(2.22)  A  cos(x+t)  +  B  sin(x+t)  +  C  cos(x-t)  +  D  sln(x-t) 
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and  for  large  positive  x   and  y  =  0  by 

A«co3(x+t)  +  B»sin(x+t)  +  C»cos(x-t)  +  D«sln(x-t)  • 

Here  the  coefficients  A^B^***  are  expressible  in 
terms  of  the  *<»#p4  ^^   accordance  with  the  asymptotic 
expressions  of  the  <p ^    (see  the  table  on  p»  32): 

B  =  -  I  Po  +  1^(^+31^2  -  I  ^2^Z 

D  =  i  Po  +  %^^-^l)'^2   -  I  ^2S 
(2,23) 

A»  =  -  K-^  ^1*^2-  1^2^  |(-^*S1)P2+  ¥2H 
C'  =  -  iV  ^1''2-  #^^3-^  ^(-^■^1)P2-  ^2^3 

Prescribing  the  incoming  wave  from  x  =  -co  is 
equivalent  to  prescribing  the  values  of  A  and  B. 
The  assumption  that  no  wave  progressing  to  the  left 
exists  at  X  =  -oo  is  equivalent  to  the  conditions 

(2.24)  C«  =  D«  =  0  • 

We  thxis  have  fovir  conditions  for  the  six  constants 

oCq»  o<'2»  *''3*P 0*^2*^3*  Consequently  the  motion  still 
depends  on  two  arbitrary  real  parameters.  Further 
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.ssinaptions  are  needed  to  determine  the  motion 
uniquely.  A  number  of  alternatives  suggest  them- 
selves here.  There  are  a  number  of  ways  in  which 
additional  physical  assxamptions  can  be  introduced 
to  arrive  at  a  uniquely  determined  motion: 

I)  No  breaking  of  waves  occurs  at  all»  Then 

This  assumption  wovild  lead  to  the  solution  found  by 
TJrsell.  For  the  present  case  of  the  vertical  barrier, 
this  may  well  be  the  solution  in  best  accord  with 
observations • 

II)  The  velocity  stays  finite  at  the  edge  of  the 
barrier: 

•^3  =P3  =  0  . 

III)  No  reflected  wave  occiors: 

C  =  D  =  0   • 

IV)  Observations  may  furnish  two  relations 
connecting  magnitude  and  phase  of  the  logarithmic 
terms  corresponding  to  breaking  with  the  magnitude 
and  phase  of  incoming  and  reflected  waves. 

A  single  condition  would  be  furnished  by  the 
assumption  that  the  velocities  are  in  phase  under 
the  barrier.  This  behavior  according  to  TJrsell  has 
been  observed  in  experiments.  It  can  easily  be  seen 
to  be  equivalent  to  proportionality  of  '^op^2  ^^^ 
tfC_,p2,  This  condition  is  automatically  satisfied 
In  cases  I)  and  II), 

In  any  case  an  ineqxiality  for  the  ^±»?'±   ^^ 
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obtained  from  the  fact  that  energy  is  either  preserved 
or  is  lost  in  the  motion.  It  can  be  shown  that  the 
rate  of  flow  of  energy  to  the  right  through  a  vertical 
plane  for  a  motion  of  the  form  (2,22)  is  proportional 
to 

A^  +  B^  -  C^  -  D^  . 

2    2 

We  can  then  take  A  +  B  as  representing  the  rate  of 

2    2 

flow  E^  of  energy  of  the  incoming  wave,  -C  -  D  as 

that  of  the  reflected  wave,  denoted  by  EL,,  and 

2     2 

A'   +  B'   as  that  of  the  transmitted  wave,  denoted  by 

Ej^  •  The  law  that  energy  cannot  be  gained  in  the 
motion  then  takes  the  form 

0  £  Ej^  +  Ey  -  E^  =  A^+  B^-  C^-  D^-  A«^-  B«^   . 


Relations  (2,23),  (2,24)  yield,  as  Is  easily 
verified, 

a2+  b^-  C^-  D^.  a»2-  b»2  =  |(SiT2-T3^S2-TrT2)(-C5P3--^P2) 

-  Tr(p62P^-P2%^  • 

This  relation  shows  that  a  change  in  energy  can  only 
occiu?,  if   ^^2*^2  ^°  ^°^  both  vanish,  l,e,  if  the 
waves  break  on  the  barrier.  If  energy  is  not  to  be 
gained,  the  constants   *^v*Pv  ™^^s*  satisfy  the  inequality 

i(SiT2-T3_S2-TrT2)e^2P3-p2'S)  ^^^^^0-^2%   ^  • 
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5,  Submerged  inrinite  barrier « 

In  this  section  we  consider  a  vertical  barrier 
extending  from  z  =  -ia  to  z  =  -loo.  We  make  the  same 
general  assiiraptlons  on  the  behavior  of  ^  or  7C  S-s 
before,  with  the  only  difference  that  X  now  can  be 
ass-umed  to  be  regtilar  at  the  origin.  The  only  singular!, 
ties  occTJir  then  at  the  edge  and  at  infinity.  This 
problem  also  has  been  solved  by  Ursell,  and  the  present 
method  will  lead  to  exactly  the  same  solutions  obtained 
by  that  author.. 

The  function  P  may  be  defined  again  by  (2,1) • 
Then  X  and  P  are  regular  in  the  region  consisting 
of  the  lower  half -plane  slit  from  -la  to  -ia\  F  is 
real  on  the  boundary  of  I  . 


The  fxinction 


z   =-~2l 


u 


maps  the  upper  w-half -plane  on  F'  ,  if   Vw  -1  denotes 
that  branch  of  the  function  that  behaves  like  w  for 
w  =  00  .  Then  Inversely 
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(3.1) 


A 


a  +z 


where  Va  +z   denotes  that  branch  that  is  imiquely 
defined  in  1""  and  reduces  to  +a  for  2=0,  Then  w 
is  an  odd  function  of  z« 


U"-/ 


U'O 


LMI 


The  condition  that 


lim  (z^  +  a'*)-J("(z)  -   0 
z  ->-ia 


implies  that 


lim  w  P  =  0  ♦ 
w->  0 


Hence  P  as  a  function  of  w  has  a  pole  of  order  ^  3 
at  w  =  0»  P  is  bovmded  everywhere  else^  and  hence 


(3.2) 


-1 


-2 


_-S 


The  function  i  -^  \  is  real  along  the  barrier, 
and  hence  as  a  fvaiction  of  w  is  real  for  real  w  with 
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,w|  <  !•  In  addition 

lim  w  i  -^  "51,=  0  • 
w->  0 

Thus  1  7^  X  9-3  ^  function  of  w  has  a  pole  of  order 
at  most  1  at  the  origin: 

a    ^^ 


i^X=g.d^w' 


Then 


P  =  (l-i^)(i4X) 


\L2n  3  ^    _00 


w 

+  (regular  terms  in  w)  • 

It  follows  that  P  expanded  in  terms  of  w  contains  no 
even  negative  powers  of  w,  and  hence  that 


cg  =  0  . 


As  a  function  of  z 


P(z)  =  Cq  +  (Cj^+Cg)  ^     -  a^C3 


VzW  y£^3 


Then 
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Vz  +a'=^ 

It  follows  that  for  large  z 

^t  +  ^y^    =  (c  +  (c,+C2))z  +  (boimded  terms), 

where  the  positive  sign  holds  for  R(z)  >  0,  the 
negative  one  for  R(z)  <  0,  As   X'  ^^^  ^^X)  are 
assvuned  to  be  bounded  for  large  z  with  Im(z)  ^  0  , 
the  relations 

°o  -^  °1  +  °3  =  °  '  °o  -  ^1  -  °3  =  ° 

must  hold.  As   "X*  +  iX   is  I'eal  for  real  z,  the 
coefficient  c,  must  be  real.  Hence,  changing  X  "by 
a  suitable  pure  imaginary  constant,  we  may  assume 
that  c.   vanishes.  Then 


x«   +  iX  = 


-  o2 


a  c 


'}[:^ 


Correspondingly 

z 


•X.(z)  =  e'^^'CaV  f  -S an+  P) 

''.CO  if?:;? 


^dn  +  p) 

.  _  ^r^^ 


z  +a  -00 


where  the  path  of  integration  Is  taken  in  T  » 


The  condition 
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Tl 

/ 

/ 

/ 
/ 

I 
\ 


v-ta 


-ta 


\. 


Im(iX»)  =  0 


for  a  point  z  on  the  barrier  then  yields 

lm(p)  =.a%lm  r  -p^zr  ^Y^ 

^-oo^^Q^+a^ 


=  -a  Cjim 


.in 


+ico  \/vf+a^ 


dY^ 


(where  the  path  rians  to  the  left  of  the  positive 
imaginary  axis) 
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-ia 


=  -a  Cglm 


2 

=  a  CgtrJ^Cla), 


Thus   Im(p)    is  uniquely  determined  in  terms   of  Cg  by 
the  formula 


Im(p)   =  a%TrJQ(ia) 


The  motion  is  thus  seen  to  depend  on  two  real  con- 
stants Cg  and  R(p)»  Those  constants  can  be  adjusted 
in  such  a  way  that  the  resulting  expression  represents 
the  effect  of  a  wave  progressing  from  x  =  -coin  the 
direction  of  increasing  x»  The  discussion  will  be 
analagous  to,  but  simpler  than  the  one  given  in  the 
last  section  for  a  finite  barrier  from  the  surface. 
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4,  Asymptotic  behavior  of  exponential  Integrals, 

In  the  following  sections  we  shall  need  some 
properties  of  the  fvmctions  H„(z)  defined  by 

(4.1)  H  (z)  =e^      /  e'  ^    F{yi)dil    (^=  0,1, ..  .,n)  . 
^         ioo 

Here  the  ^^     are  defined  by  (1.13).  F  is  regular 
in  P '" ,  real  along  the  real  axis,  and  has  a  zero  of 
order  (n  +  1)  at  infinity,  as  indicated  by  (1.25), 


The  integrals  in  (4.1)  converge  because  of  (1.14). 
H„(z)  is  uniquely  defined  in  the  region  "P  '"'  slit  along 
the  negative  real  axis,  if  the  path  of  integration 
starts  out  along  the  lower  border  of  that  slit.  They 
satisfy 


(4.2)  H»(z)  -r^  Hjz)  =  P(z) 


The  values  of  H  (z)  on  opposite  borders  of  the 
negative  real  axis  are  connected  by  the  relation 
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(4.5)    H„(z  +  oi)  =  H  (z  -  lo)  +  C^  e  ^ 


valid  for  negative  real  z,  where  C^  denotes  the 
constant 

(4.4)  C^  =  dr  e"  ^     P(t^)dY( 

(the  integral  being  extended  over  a  circle  of  radius 
>  a  about  the  origin).  Thus  the  function 

(4.5)  H^(z)  -  ^  C^  e  ^   log  z 

is  xini valued  in  T*  ""*. 

The  H^  satisfy  the  symmetry  condition 

(4.6)  H^(z)  =  \      iz)   +  G^e  ^ 


Proof  of  (4.6); 

Prom   (1.16)    it  follows   that  for  negative  real  z 


jz 


e 


'-00  "-00 


{^)di{=         e     ^/"p(v^)dyj 

•/-no 


if  the  path  of  integration  is  taken  along  the  negative 
real  axis.  Thus 


v^^  ■  w^ 


for  z  on  the  lower  border  of  the  negative  real  axis. 
Continuing  analytically  across  the  negative  real  axis 
the  left  hand  side  goes  over  into 
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H^{z)  -  C^e  ^ 

according  to  (4,3),  the  right  hand  side  goes  over 
into  the  analytic  fiinction 


As  both  those  functions  are  analytic  and  univalued 
in  V  *  slit  along  the  negative  real  axis,  they  must 
be  identical  throughout  that  region. 

Replacing  jx   by  n-^  in  (4,6)  we  obtain 

Comparing  vrith  (4.6)  and  using  (1,16)  it  follows  that 


(4.7)  C^  =  -<^n-/^. 


Lemma  I; 

If  Ca^  =  0,  the  function  BLjz)  is  regiilar  at 
infinity  and  has  a  zero  of  order  n  +  1  there, 

^roof : 

C^  =  0  implies  that  Hu.  is  univalued  in  K'"'. 
Then  H  .  permits  a  Laurent  expansion  of  the  form 


H  (z)  =  ^_  YvZ 

'  k=-'Co 


for  |z|  >  a.   Then  according  to  (4,2)  and  (1,25)  the 
Yv  satisfy  the  recursion  formula 

(k-l)Yjj..i  +  ^/^  Yjj.  =  -«^  (-00  <  k  <  +00  )  . 
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As   <^j^  vanishes  for  k  ^  n,  it  follows  that 
Yi  »  Yg  =  ...  =  Yj,  =  0 

and  that 

-k 


Yk 


=  -^  Yo  :f  or  k  <  0 


Hence 


r  k=n+l 


Here  Y  must  however  vanish,  as  obviotisly  by 
'  o 

definition 


lim  H^(z)  =  0  . 


Lemma  II: 


For  I z I  large 

O(-l^)  for  R(z)  <  0,  I^(z)  <  0 


(4,8)   H^(z)  = 


C  e"   +  0( — \r)   otherwise 
z 


Proof: 


We  have 

Z  2 

(4.9)  EL(z)  =  e"^^  \     e^^l  F(y\,)dYV=  -iP(z)+le-^^  f  e^^P»(^)dr( 
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According 


to  (1.25)  there  Is  an  M  such  that 


|p«(z) 


!!n+S 


For  -  2  <  9  <  0  we  can  choose  the  path  of  integration 

so  as  to  rvin  from  +ioo  to  ri  along  the  positive 

i© 
imaginary  axis,  and  then  from  ri  to  re   counterclock- 
wise along  the  circle  |y|_|  =  r.  Along  this  path 
|e^|  <  le"^  ^|«  Consequently 

|H^(z)|  <  |F(z)l  +f  -4^^A 

n 

<|Hz,|.M(_l-^.^,=0(-i^, 


Which  proves  the  first  part  of  (4,8),   If  0  <  ©  <  -iji  , 
we  split  the  path  into  two  parts.  The  first  part 
shall  run  from  +ioo  arotind  the  circle  \^\    =  r  back  to 
+ioo ;   the  second  part  from  +ioo  along  the  positive 


imaginary  axis  to  ri,  and  then  along  the  circle 
|yi|  =  r  along  the  shortest  path  from  +ri  to  re   . 
The  integral  over  the  first  part  has  the  value  C 
On  the  second  part  | e  ^|  <  | 
the  same  estimate  as  before. 


+iZi  ^ 

e   I >  so  that  we  have 
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5,  Derivation  of  the  complete  system  of  conditions 
on  the  general  solutions  for  the  inclined 
barrier. 

In  section  I  it  was  seen  that  the  most  general 
function  P  satisfying  the  boimdary  conditions  for  P 
and  having  the  proper  singularities  is  represented  by 
(1,25),  and  thus  depends  on  5n  -  2  arbitrary  constants 
Cv«  The  boundary  conditions  for  P  are  however  not 
equivalent  to  the  original  conditions  (1,6),  (1.7)  for 

")(  (z),  which  have  to  be  satisfied.   The  conditions 
on  the  c,^  that  guarantee  that  (1,6),  (1,7)  hold,  will 
be  derived  now. 

The  rational  function  1/l(D)  admits  a  partial 
fraction  representation  of  the  form 


1    _  '^  "^k  ,  ^  ^k 

L(D)        k=l  ^~^k       k=l  D^ 

• 

_          1                         7-1 

^        l(d) 

LMtk)     '      ^       ^^-^^^ 

LdD--^ 

(5.1) 


Here 


(5.2) 


More  explicit  expressions  for  the  CT      and   T, 
can  be  obtained.  Thus  e.g. 

(5.3)  <r^  =  1(1^^)^-1  i  . 


Por  what  follov/s  only  certain  symmetry  properties  of 
the  0*^  and  Z^  corresponding  to  (1.16)  and  (1.17) 
are  of  importance. 

We  have  from  (1.19) 
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where  Q'  is  a  polynomial  with,  real  coefficients • 
Hence,  using  (1.17) 


(5.4)  e<y,   =   €3: 


k  "  ^  n-k+1  ' 


Similarly  it  follows  from  (5.2),  (1.19)  that 


(5.5)  Im(e^Tv)  =  0  • 


We  have  from  (1.18) 

!:=1  ^    ^         k=l 


As  l/P(D)  is  a  rational  function  with  real  coefficients 
and  vanishing  at  infinity,  it  follows  that 

n 


(5.6)  ^  ^k  +  ^1=0 

k=l  ^     ^ 


(5.7)  {i;,y+l)<T^  =  i^^^^+^)<y^^^     for  k  =  l,...,n-l 


* 


Relations  (5.3),  (5.7)  for  the  (T.    are  similar  to 

those  found  for  the  coefficients  c.  in  the  paper 

by  J.  J.  Stoker,  AMG-NYU  153,  p.  95.   They  could 
be  used  for  explicit  calculation  of  the  CT.  • 
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,5,8)  ^('^k+l  "*  ^''^k^  ^  ^  ^°^     "^  ^   l,...,n-l 


(5.9)  Im(ir^)  =  0  . 


In  accordance  with  the  partial  fraction 
resolution  (5,1)  v;e  can  write  the  general  solution 
of  the  differential  equation  (1,8)  for  X  (z)  in  the 
form 

X(z)  =  L(D)-^  P 


k=l  ^    ^       k=l  ^ 


^1-     Cj^z   n=i    , 
k=l  ^       k=0  ^ 


Here  the  a,  and  b^  are  arbitrary  constants, 
(D-i^v)"'^  denotes  a  special  fiinction  with  the  property 
that  the  operator  i^-Cir)    applied  to  it  gives  P,  and 
D"T'  a  special  function  with  P  as  k-th  derivative.. 
We  are  free  to  choose 

(5.10)   (D-^i^)  -^P  =  e  ^  1   e"  ^  P(^)dY^  =  Hj^(z) 

^-00 


/z 
(z-Y|)^-l  P(Y()dy|  =  h^(z)  . 

-00 


Because  of  (1.19)  and  R(4k)=  ^*   ^11  integrals  converge 
absolutely. 
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As  H^  (2)  and  h^iz)    tend  tov;ards  0  for 

z  -S"  -{x>,  and  as  "X.Cz)  is  bounded  at  Infinity,  it 
follows  that  the  real  part  of 

n_     Z,,z       n-1     ^ 

k=I  ^       k=0  ^ 

must  be  bounded  for  z  ->-co •  As  all  the  ^-  with 

the  exception  of  ^   have  negative  real  parts  and 

are  distinct,  it  follows  that  all  a^^  and  bj^  with  the 

exception  of  a  ,b  .b,  must  vanish.  Then  for  large 
n  O''  1 

negative  real  z 


(D  +  i)  X(2)  -  (bi  +  ib„)  +  ib.Z 


and  hence  because  of  (1,24) 


^o  =  ^1  =  °  • 


Consequently,  writing  oC  for  a_, 

(5.12)  X(z)  =ode-^^  +X<J-kHk(z)  +^  "^k^^^) 

K=l         k=l 

Y/e  have 
(5.15)  ^  ^k  =  ^k^k  ^  ^ 

(5,14) 

P    for  k  =  1 

and  hence,  using  (5,6), 
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f  =  (D+1)  X  =  ^  ^k^^k^^^^k^^^ 
(5.15) 


According  to  (l.ii4)  tiie  function 
f  =  (D  +  i)  X 
is  regular  at  infinity.  On  the  other  hand  (4.3)  shows 


that  Hw  (z)  changes  by 


Cj^e 


^k^ 


along  a  closed  path  surrounding  oo  .  The  functions 
h.  (z)  are  regular  at  infinity;  indeed  from  (1,25) 
they  have  an  expansion  of  the  form 


(5.16)   h,  (z)  =  d"^  =  2 


-k„   -^—     ""k^ 


k-r 


r=gi+l 


(k-r)(k-r+l)...(l.r)   * 


Consequently 


g  ^k^^k-^^)  ^k-    =°  • 


It  follows  that 


(5.17)   Cj^  =  d)  e  "^  P(T[)dY|^  =0  for  k  =  l,,,.,n-l 


where  the  integrals  are  extended  over  a  path  with 
|y||  >  a  . 
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Condltions  (5.17)  are  already  sufficient  in 
order  that  X  (2)  ^s  the  required  behavior  at 
inl'inlty  and  along  the  real  axis.   We  have  Indeed 
from  (5,12),  (5,16),  (4,8)  for  large  |z|  with 
R(z)  <  0 

(5.18)  X(z)  =  '^e-^''  +  0(|) 

and  from  (5.17),  (4,8)  for  large  |z|  with.  R(z)  >  0 

(5.19)  X(z)  =  (oC+  <r^C^)e-^'^  +  0(i) 

with  corresponding  formulae  for   ^'(z).  Thais   ')(_*(z) 
is  bounded  at  infinity  in  T  ,  and 


lim   X*(2:)  =  0 
y->-oo 


It  follows  from  (4.6),  (5.17)  that  for  real  z  ;^  0 
hJTzT  =  Hj^.k(z)  i'oi'  k  =  l,...,n-l  , 

Besides  J1t,(z)  Is  real,  as  the  coefficients  oc,    in 
(5,16)  are  real.  Then  (5,7),  (5,8),  (5,9)  show  that 
all  terras  in  the  representation  (5,15)  of  f  are  either 
real  or  conjugate  imaginary.  Consequently  f  =  (D+i) X 
Is  real  for  all  positive  and  negative  real  z  f^  0» 

It  remains  to  derive  the  restrictions  on  P 
introduced  by  the  barrier  condition  (1.11).  For  this 
piirpose  we  introduce  the  abbreviation 


(5.20)    i^(z)=^i^X  ^or        /x,=  l,...,n 
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Here  L(D)/(D-^)  Is  of  course  a  polynomial  in  D. 
We  have  from  (1.12) 

(5.21)       (D-^^)  r^(z)  =L(D)X  =  P(z)  . 

Hence,  as  H^  is  a  solution  of  the  same  differential 
equation, 

y  (z)  =  H^  (z)  +  d^e  ^ 

with  a  certain  constant  d^,  •  The  constants  di^^ 
can  be  determined  from  the  behavior  at  infinity. 
X  and  its  derivatives  are  bounded  for  real  z  with 
large  absolute  values;   so  is  H^(z)  on  the  basis 
of  (5.17).  Hence 

dyLA.  =  0  for   /^=  l,...,n-l  . 

For  large  negative  real  z  all  Hi^,hT^(z)  and  their 
derivatives  tend  to  zero,  and  hence  from  (5.12) 


viz)   ^   ^  oce-^^  =  L'(-i)  06  e-^^ 


^^      ^^ 


=  0-^  cxe 


Consequently 


(5.22)        r^(z)  =H^(z)  +^    8^  e-^^ 


We  now  form  the  expression 
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c-^{^      ^  ^  \    -    c  "l^L(D)  ^  _L(D)   V  -y 


As  Q  has  real  coaTficients  and  (1.13)  holds,  the 
differential  operator  applied  here  to  (6DX)  is  a 
polynomial  in   6D  with  real  coefficients.  Hence, 
when  (1,11)  holds,  it  follows  that 

(5.23)  6-^(r^  +  v^^^^)   =  u 

is  real  along  the  barrier.  Applying  the  operator 
(6D)  to  u  we  again  get  a  real  expression 

(see  (5.21))  and  hence,  as  P  is  real,  the  expression 

(5.24)  2:^r^+  ^n.^A.^1 


=  V 


is  real  along  the  barrier.  If   ll  j^  n-/u+l  we  get 
from  (5.23),  (5,24) 


n-^1      ^     ^^n-/*H  V^ 
and  thus  with  the  help  of  (1.17) 
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i  e  -^p^  =  -i  e  -\.^+i 


or 


The  same  relation  follows  from  (5,23)  alone,  if 
/*.  =  n-u.+l.  Substituting  the  expression  (5.22) 
for  r^  ,  we  see  that  along  the  barrier 

(5.25)  H^(z)  =  6^  *in._^i(^)  m|:S^-^  fi^je^*^" 

n 


for  /x-  =  l,,.,,n 

Similarly,   putting 

we  find 


^      X=  St,(z)    , 


D^  Sjj  =  P   • 


Hence  s,^  and  h-  differ  at  most  by  a  polynomial.  As 
both  3-^   and  h^  vanish  at  infinity,  they  must  be 
identical:   Then 

e-\(z)  =^7C  =^  (£DX)  . 

As  Q(eD)/(eD)   for  k  <  n  is  a  polynomial  in  (6D)  with 
real  coefficients,  it  follows  from  (1.11)  that  along  the 
barrier 
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(5.28)     Ira(  6  "-^1,  (z))  =0  for  k  =  l,...,n-l 


^k^ 


Conditions  (5,26),  (5,28)  which  have  been 
derived  fi->om  (1,11)  are  actually  equivalent  to 
(1,11),  as  can  easily  be  verified  from  the  expression 
(5,12)  for  X   ,  if  use  is  made  of  (5,4),  (5,3),  (1,17), 

The  complete  set  of  conditions  on  P  consists 
then  of  (5,17)  and  the  relations  (5,26),  (5,28),  which 
have  to  be  satisfied  for  every  z  on  the  barrier.  We 
shall  replace  the  identities  (5,26),  (5,28)  by  a 
finite  set  of  conditions  in  the  next  section. 


.../^ 
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3.  Complete  finite  systeia  of  conditions  for  the 
solution. 

In  order  to  permit  a  reduction  of  the  conditions 
for  P  to  explicit  equations  for  the  coefficients  Cj^, 
it  is  convenient  to  study  the  conditions  in  terms  of 
the  variable  w  Introduced  in  (1,31).  We  denote  by 
Tt.  for  k  =  0,,,.,4n-l  the  sector  defined  by 


k        ^  (k+1) 
^  <  arg  w  <  i.^ 


We  define  the  points  A^  and  B^  of  the  w-plane 
respectively  by 

(6.1)      A^  =  e  -^-2^  ,  B^  =  (2n-l)-  ^  S'^^ 


Then 


Ao  =  ^"^   '   B^  =  /^ 


are  the  points  on  the  boundary  of  T  corresponding 
respectively  to  the  edge  of  the  barrier  and  to  the 
point  at  infinity  of  the  z-plane.  The  fxmction  G(w) 
is  real  on  the  boundary  of  T  and  hence  on  all  rays 
OA  00  and  OB^oo  ,   It  satisfies  the  relation 

(6.2)  G(  e^w)  =  G(w)  . 


G  has  poles  of  order  ^  4n-l  at  the  A  ,  zeros  of  order 
(n  +  1)  at  the  B  and  poles  at  0  and  oo  .  The  rational 
fvinction  K(w)  has  poles  of  order  1  at  the  B  • 
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It  is  then  obvious  that  the  functions  Hj^  and 
h-  as  f\inctions  of  w,  which  are  integrals  of  uni- 
valued  fvinctions,  can  be  continued  into  the  whole 
w-plane,  with  possible  branchpoints  at  the  points 
A  ,B  ,0,00  •  The  same  holds  then  for  the  ftinction  X  • 
It  will  first  be  shown  that  as  a  consequence  of  the 
barrier  condition  no  branchpoints  actvially  occur  at 
the  Ap. 

We  observe  that  the  function   fX^Cz)  is  real 
along  the  barrier,  and  satisfies 


(6.3) 


lira  (z-ae)  e"X»(z)  =  0  , 
z  ->•  ae. 
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""lie  transformation 


(6.3a) 


z  =  6  (a-^'^') 


maps  the  interval  -\/a  <  C   <  +  sf^     o^  tlie  real 
K,   -axis  on  the  two  borders  of  the  barrier  in  the 
2-plane,  Then  6X  '  as  a  function  of  C    is  I'eal 
along  the  real  ^  -a:cis  in  a  neighborhood  of  the 
origin,  and  hence  lonivalued  in  a  full  neighborhood 
of  the  origin. 


Besides  because  of  (6,3) 


lim  ^^     tX'  =  0 


so  that   e  %  '  ^s  a  pole  of  order  at  most  1  at 
C  =  0,  Consequently  an  expansion  of  the  form 


(6,4)   £X»  =^  ^^^=^   aj^(a-&-^z)^/2 


k=-l 


k=-l 


is  valid  near  the  edge.  As  P  is  formed  by  applying 
the  (2n  -  l)st  order  differential  operator  Q(6D)  to 
gX'#  ^^   follows  that  P  permits  an  expansion  of  the 
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(6.5)     P(z)  =  ^~  b.(a-e-^)^/2  =;^Z  V^  . 


It  is  obvious  that  no  negative  Integral  powers  of 
(a-e"~z)  can  arise  by  differentiation  of  (6.4), 
Hence  we  have  that  In  the  escpanslon  (6,5)  of  P  the 
coefficients  with  even  negative  Index  must  vanish: 

(6.6)  b^2  =b_4  =  ...  =b2^^  =  0  . 

Let  now  U(z)  denote  an  arbitrary  fxinction  of  z, 
which  Is  regular  at  z  =  at.  Then  relations  (6.6) 
imply  that 

(6.7)  U(z)  P(z)  ^ 


as  a  function  of  ^  has  residue  0  at  C    -  0»     Indeed 
the  expression  (6.7)  reduces  to 

52—    k+i 


-2  eTKea-e^^**)  >   b.  c 

k=l-4n  ^ 

and  no  odd  negative  powers  of  ^  can  occur  in  the 
expansion  of  this  expression  with  respect  to  C  , 
if  (6,6)  holds.  On  the  other  hand  relations  (6.6) 
can  be  seen  to  be  equivalent  to  the  conditions  that 
the  residue  of 

(6.8)      (z-ae)^  P(z)  ^   (k  =  0,l,...,2n-2) 
considered  as  a  fvoiction  of  J  vanishes  for  ^  =  0. 
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CoEiparing  (6.3a)  with  (1,31)  we  have 

^2  -  o   (2n-l)w^-l-2n£"V^"^ 

r    -  t ^ — • 

^  (2n-l)w'^^-l 


Here  the  numerator  on  the  right  hand  side  is  a 
polynomial  in  w  which  has  a  zero  of  order  2  for 
w  =  6"  ,  It  follows  that  in  the  neighborhood  of 
w=  t"     or  C=0  the  variables  ^  and  w  are  imi- 
valued  fimctions  of  each  other.  Consequently 


f 


U(z)  P(z)  ^  d^ 


extended  in  the  t,  -plane  around  a  path  aroimd  the 
origin  goes  over  into 


f 


n(K(w))  G(w)  Kt(w)  dw 


over  a  closed  path  about  the  point  w  =  6  ,  Thus 
the  vanishing  of  the  residue  of  (6,7)  at  ^  =  0  is 
equivalent  to  the  vanishing  of  the  residue  of 

(6.9)  U(K(w))  G(w)  K«(w) 

at  w  =  &"  .  The  vanishing  of  the  residue  of  (6.9) 
for  any  regular  fimction  U  is  equivalent  to  the 
vanishing  of  the  residue  of  the  f;mctions 

(6.10)  (K(w)-ae)^(w)Kt(w)   for  k  =  0,l,...,2n-2  . 


We  have  then  the  condition  that  the  coefficient 
of  (w-e*  )"  in  the  expansion  of  (6,10)  in  powers  of 
(w-e**  )  vanishes.  Using  the  expression  (1,29)  for 
G(w)  this  obviously  yields  (2n-l)  linear  homogeneous 
equations  for  the  Cj^  with  coefficients  that  are 
rational  numbers,  independent  of  a. 
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2,v 

Replacing  w  in  (6.9)  by  few  results  in  the 


expression 


6"'^U(e"^^K(w))G(w)K«(w)  , 


which  is  again  of  the  same  form.   It  follows  that 

for  a  fimction  U(z)  which  is  regular  at  all  points 

— ^r 
2=6"  the  expression  (6,9)  has  residue  0  at  all 

points  Ap. 

By  the   substitution   (1.31)    the  fvinctions  R,    and 

h,^  defined  by   (5.10),    (5,11)    go  over  respectively 

into  functions 

4iJC(w)    p^      .?T^(V) 
(6,11)   Pj^(w)    =  Hj^(K(w))    =  e  ^  /      e     ^        G(v)Kt(v)dv 


(6,12)  Pj^(w)    =h^(K(w))    =jErrTr   /   ^^^""^    -K(v)]^"^GK»dv 

^o 

where  the  path  of  integration  starting  with  Bq  first 
rxms  along  the  real  axis  to  the  left. 
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Pt^(w)  and  Pj^(w)  are  then  uniquely  defined  fimctions 
in  T  and  can  be  continued  outside  that  sector  by 
means  of  formulae  (6.11),  (6.12)  with  values  depending 
possibly  on  the  choice  of  the  path  of  integration. 
The  functions  of  Yl 

e  ^    and   (z-tj)^"^ 

are  regular  for  all  V^  ♦  Talcing  them  as  U  in  (6,9) 

it  is  seen  that  the  Integrands  in  (6.11),  (6,12)  have 

residues  0  at  all  A^,  and  hence  that  the  Pi,(w)  and 

r'  Ic 

Pj^(w)  have  no  branchpoints  at  the  A   (whatever  path 
may  have  been  chosen  in  reacliing  those  points). 


To  define  T^   and  Pj^  in  T^_^-,  in  agreement  with 

the  convention  made  already  in  defining  Hj^  and  h.  in 

"T*  ''»   we  have  to  take  the  path  of  integration  so  as 

to  cross  the  real  w-axis  in  the  interval  between  B^ 

o 

and  +00  , 

V/e  first  discuss  the  behavior  of  the  functions 
Pjj(w) ,   The  integrand  in  (6,12)  has  a  zero  of  order 
n-k  at  V  =  Bq,  and  hence  is  finite  at  B  ,  As  G(v)  and 
K(v)  are  real  for  real  v,  it  follows  that  Pij.(w)  is  real 
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for  real  w.  Besides,  because  of  (6,2)  and 
(6.12a)  K(e^w)  =  £"nC(w)  ,  K»  (£^w)  =  e"\t(w) 
we  have 

(6.13)  Pj^(e^w  =  6-2^j^(w)-  6-^^J   (K(w)-K(v))^"^G(v)Kt(v)dv 

B_ 


As  K(w)  is  rational,  it  follows  that  Ptj.(w)  has  no 
branchpoints  at  any  of  the  points  B  • 
We  observe  that  at  infinity 


(6.14)      G(w)  =  C5n«2''   ■*■  (^o'^^'^^^i  terras) 


Similarly  at  infinity 


K(w)"^  =  (^2i)2n  ^2n  ^  (bounded  terms)  . 


Hence,  denoting  by  f   the  constant 


Y  -  ^^n=l)    °5n-2  ' 


we  have  in 


(6.15)        g(w)  =  G(w)  -  r  K"^(w) 


a  f "unction,  which  is  bounded  at  infinity.  As 
K(v)  =  0(i)  K«(v)  =  O(^) 


we  have  for  w  =  t  t       with  t  >  1  , 
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/CD 
[K(w)-K(v)]^-^g(v)K»  (v)dv 


(6.16) 


+   (kll);   /  [K(w)-K(v)]^'V^(v)Kt(v)dv 
B_ 


■'•T^ 


00 


e-H 


,k-l 


[K(w)-K(v)]'^"-^g(v)K«(v)dv  , 


Here  the  second  integral  has  the  value 
pK(w) 

-00 


(l-Sn):.:(k-2n)-  -  (i-J>n)..:(k-§n)   • 
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The  third  integral  in  (6,16)  is  of  the  form 
t    (real  number),  as   fe^T^Cw),   t^TCCv),  g(v), 
e<iv,  e  K»(v)  are  all  real.  Now  according  to 
condition  (5,28)  E"  ri.  (z)  must  be  real  along  the 
barrier,  i,e,   6"  i>Tj.(w)  must  be  real  along  the  ray 
0,  6"  ,00 ,   It  follows  that,  as  it  is  already  known 

that  Pv(w)  has  no  branchpoint  at  t"   ,  it  is  sxofficient 

-k  -1 

to  have   fc  i>T,(w)  real  along  the  segment   6   »oo , 

— k 
The  reality  of  Pv.(w)  t"     on  that  segment  is  according 

to  (6,16)  equivalent  with  the  reality  of 
(6.17)    £"^  /   [K(w)-K(v)]^-^  g(v)K»(v)dv 


Bo 

for  w  =  e   t  with  t  >  1,  The  e35)ression  (6,17)  is  a 
polynomial  in  the  real  variable   6~"T^(w)  =  K(t), 
All  coefficients  of  that  polynomial  have  to  be  real, 
l,e, 

e"*^"^  /   K(v)^~^  g(v)Kt(v)dv   (r  =  0,,,,,k-l) 

Bo 

must  be  real.  As  the  integrand  is  real  along  the  path 

of  integration,  which  falls  into  the  real  v-axis,  and 

—l—r 
£     is  not  real,  this  amovints  to 


,00 


<6.18,    /  K(v)'-s(v)K.(v,d.  =  0  for  r  =  0 n-S 


(k  being  any  of  the  niimbers  1,  ,,,,n).  Relations  (6,18), 
which  are  seen  to  be  equivalent  to  (5,28),  constitute 
(n-1)  homogeneous  eqxiations  for  the  c- ,  with  coefficients, 
which  are  real,  are  independent  of  a,  and  can  be 
evaluated  explicitly  as  definite  integrals  of  rational 
functions. 
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Incldentally  It  is  seen  from  (6,16)  that 
p,  (w)  has  no  branchpoint  at  w  =  oo  ,  as  the  first 
integral  has  been  seen  to  vanish,  the  second  one 
is  a  rational  function  from  its  explicit  expression, 
and  the  integrand  of  the  third  one  is  0{-^)    at 
V  =  00  ♦  The  function  Pt^(w)  having  no  brXnchpoints 
a-c  the  points  A.  #B  ,oo,  cannot  then  have  a  branch- 
point at  the  only  other  possible  place  w  =  0,  As 
p,  (w)  Is  formed  by  integration  of  rational  functions, 
it  follows  that  Pj^Cw)  is  a  rational  function  of  w, 
with  its  only  possible  poles  at  the  A^  and  the  origin, 


We  pass  on  to  the  consideration  of  ^he  functions 
P,  (w),  which  have  been  defined  uniquely  in  the  region 
formed  by  the  sectors  T^  and  T^^t  ,  slit  from  0  -co 
B  along  Che  real  axis.  This  region  is  the  image  of 

V ""'  slit  along  the  negative  real  z-axis  -under  the 
mapping  z  =  K(w).  Conjugate  imaginary  points  go  over 
into  conjugate  imaginary  points  under  that  mapping. 
Then  relations  (4,6)  yield 
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(6.19)   P^(w) 


=  W^^  ^  ^^^' 


^u.K(w) 


for  yU/  =  0, , 


»  •$! 


These  are  not  conditions  for  the  c-^   since  they  are 
satisfied  Identically, 

As  a  consequence  of  the  conditions  already 
obtained  from  (6,10),  the  Pt-(w)  have  no  branchpoints 
at  the  A  .   Then  Eu(w)  is  iimnediately  defined  iiniquely 
in  the  region  consisting  of  the  four  sectors  T^  o^ 
^4n-l'^o''^l'  sl^*  from  0  to  B  .  By  analytic  continiiation 

relation  (6,19)  persists  in  this  larger  region,  Apply- 

-2 

Ing  the  substitution  v  =  t     u  in  (6,11)  we  have,  using 

(6.2),  (6,12a) 


Pk(w) 


^^C^.-,K(e^w)^2 


e 


-^k-lK(^) 


G(u)K»(u)du 


^2n-l 


where  the  path  of  integration  lies  in  the  region 

formed  by  '^4^^4,>'^4^^5»'^4:ri~2''S^^l  ^^^^   ^^°^S  OBg^^^^' 
and  rims  first  along  the  upper  side  of  the  segment 

OBg  n .   If  now  k  >  1,  and  if  w  lies  In  T  +  T,,  the 

path  of  integration  can  be  broken  up  into  a  part 

running  from  Bg^^-i  to  B  and  a  part  rxmning  from  B 

to  6  w. 
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Tliere  results   the  formula 


(6,20) 


Pjj.(w)   =  e^p 


j^^^(e^w)+  t^e  ^  \         e     ^  -^         G(u)K» 

^2n-l 


du 


2  2  ^k^^^^)     A  -^k^(^) 

=  e'=^Pj^_^(e'^w)    +  e  ^  /       8     ^        G(v)K»(T)dv 


(k  =  l,..,,n)    . 


( 


Here  the  path  of  integration  in  the  last  integral 
starts  along  the  upper  side  of  the  segment  OB  and 
ends  on  the  upper  side  of  the  segment  OB, .  We  have 
from  (6,19) 


Vl(^^^)  =  fn.k+l(^)  +  Cj^-l  ^ 


Cj^K(w) 


so  that  (6,20)  may  also  be  written 


-71- 

"B 


(6.21)  Pj^(w)  =  t"^Pn-k+l^^^^"^  ®      ^^k-l"^J   ®       GKtdv) 


o 


for  w  in  T  +  T,  and  k  =  l,.,.,n.  Replacing  k  by 
n-k+1  and  w  by  £^w  (which  again  lies  in  T^  +  T^) 
we  get  from  (6.16) 

W.l(^"^)  =''^^  ^   e^-^^^\c,.,./'e-^-^<K.dv]  . 


Comparing  with  (6.21)  and  using  (4,7)  it  follows  that 


(6.22)  f    V^^-^-^1^  ""  G(v)K»(v)dv  =  Cj^  -  e,''\^^ 


Bo 


/ 


Bn 
B, 


-  e"*^  /   e'  ^    G(v)K»(v)dv 


for  k  =  1, . . .,n. 


We  are  now  in  a  position  to  formulate  conditions 
(5.25)  in  a  simple  fashion,  z  on  the  barrier  corresponds 
to  w  =  f^t   with  real  positive  t.  For  such  w  the 
condition  is  that 

(6.25)  Pi^(w)  =  t\^^^^{y)   +  (^  8^^  e^  -  ^  Sg)e  ^    . 

n 


For  w  =  e**"^t  we  have  w  =  6  w  .  Thus  using  the 
identities  (6.21),  relations  (6.23)  are  seen  to  be 
equivalent  to 


-72- 


(6.24)  /   e-  ^    G(v)Kt  (v)dv  =  -Cj,  n-^-  6^  e'^-  f-  5^ 

k  =  1, . . . ,n  • 


We  have  in  (6,24)  a  system  of  n  linear  equations 
with  complex  coefficients  for  the  real  numbers  Cj^« 
In  (5.17) 

(6.25)  ^k  ~  °  ^°^  ■^  ~  1»  •••»"-! 

we  have  another  such  system  of  n-1  eqimtions.   In  view 
of  the  identical  relations  (4.7)  and  (6.22)  these 
2n-l  equations  with  complex  coefficients  can  be  seen 
to  be  equivalent  to  the  (2n-l)  equations  with  real 
coefficients: 

(6.26)  R(&"-'-Ct-)  =  0  for  k  =  l,,..,n-l 


f^l   .^iJC(v)  -   1  2  oc   n> 

(6.27)  R(  y   e"  ^    G(v)K»  (v)dv  +  Cy^^^-  ^  S^  6^^+  ^  ^  J)  =0 

for  k  =  l,,.,,n  . 
Indeed  (6,26)  implies  that 

t"   Cj^  =  ivjj.  for  k  =  l,,..,n-l 


with  real  Vj^,  Then  also 


^"^^n-k  =  i^n-k 
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and  hence,  using  (4,7) 


^"'"Cn-k  =  -^Cj,  =  -i^n.k  • 


Thus 


fe  ^k  =  ^n-k 


2 

and  hence,  as  £   Is  not  real  for  n  >  1, 


^k  =  ^n-k  =  0  . 


Similarly,  using  (6,22)  relations  (6,27)  already  imply 
(6.24), 

V/e  have  found  now  that  the  boiondary  conditions 
for  A.  are  equivalent  to  the 

(2n-l)  +  (n-1)  +  (n-1)  +  n  =  5n-3 

linear  equations  with  real  coefficients  for  the  c-^ 
respectively  obtained  from  (6.10),  (6,18),  (6,26), 
and  (6.27). 

These  relations  can  be  written  explicitly  as 
linear  equations  for  the  c,^,  if  we  introduce  the 
rational  fvinctions 

r«2n   1  iH+l 
t"  -  ^1=1^         2nr 
^r(^)  =  K'(w)  ^2n(Sn-l)(3^^^2n)4^>l  "^ 

for  r  =  0,1,. , ,,5n-2.  Then  (6,10)  can  be  written 
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5X1^     r  V 


m-^ 


(K(w)-  ae)""  4p(^)«i^  =  0 

for  k  =  0,l,...,2n-2 

where  the  integral  is  taken  over  a  closed  path  about 
oint  S"   • 
Relations  (6.18)  can  be  v/ritten 


the  point  S"   , 


>"  c^  f  K^(w)  4_(w)dw  + 
k=0   ^  J,.. 


k=0     v^ 


-5n-2  /"K^(-)t4  5n-2(-)  "  ^  (2n-l)K(w)  )'"^^"  =  ° 


for  k  =  0,1,.,., n-2  . 


Relations  (6,27)  become 
5  n-2 


Cy  R(&"-^d)  e"  ^    4p(w)dw)  =  0 


r=0 

for  k  =  l.....-n-l 

where  the  integral  is  extended  over  a  closed  path 
abouc  the  point  B^  =  /^.  Finally  (6,27)  becomes 

5n-2        n  ^     V^(w)  ^    -,2 

^ZICpRt      e-^    4(w)dw.^sj£2 
^=0       J^  Oi, 


^  S2  +  i  e-^-^     4p(w)dw]  =  0 


for  k  =  1, ,.,,n 
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where  the  path  for  the  first  integral  is  represente.l 
by  Pig.  24  ,   and  for  the  second  integral  is  a  closed 
path  about  yU/  • 

We  have  then  (5n-3)  linear  equations  for  the 
(5n-l)  quantities  Co'***'°5n-2  ^  terms  of  the  coniplex 
nximber  oC  ,  which  describes  the  behavior  of  X  (z) 
at  X  =  -00,  We  can  still  prescribe  the  strength  of 
the  energy  sinks  at  the  two  corners  of  I  at  the 
origin.   This  amounts  to  prescribing  the  values  of 

c^  and  Cc  «»  the  vanishing  of  which  indicates  the 
o      on— 2  °  _. 

absence  of  logarithmic  singularities  of  "X  at  those 
comers.  Assuming  c  and  Cgj-^o  ^o  be  prescribed,  we 
then  have  5n-3  equations  of  the  general  form 

(6.28)   %(Ci*...»C5^.3)  ^^ii<^o*Hn^2'^''^^ 

for  i  =  l,,..,5n-3, 

where  the  L^  and  Mj  are  homogeneous  linear  functions 
of  their  respective  arguments,  and  where  we  put 

p  =  RM      »   Y  =  Im(e6)  . 

We  shall  not  prove  here  that  the  eqiiations  (6,28)  are 
linearly  independent  and  hence  have  a  unique  solution, 
but  we  shall  prove  the  corresponding  property  for  a 
derived  system  of  equations,  that  takes  into  account  the 
time  dependence  of  the  function  $  • 

Under  the  assumptions  made  in  section  1  the 
velocity  potential  of  the  motion  of  the  liquid  is  of 
the  form 

$(x,y,t)  =  (p»(x,y)cos  t  +  9"(x,y)sin  t  , 
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v/here 

and  the  fiinction  7( »  and  X"  separately  satisfy 
all  the  conditions  Imposed  above  on  X^^).   (Here 
the  symbols  "^^   and  ")(,"  ^^®  used  to  denote  inde- 
pendent functions  and  are  not  to  be  confused  with 

derivatives  of  a  fimction  X  •)   If  we  denote  by 

t   t   t       II   II   ii 
Ci^>P  ff     and  c-i^,p  ,Y  respectively  the  coefficients 

Ct^  entering  the  expressions  for   X'  ®^^  X"'  ^® 

have  the  lOn-6  equations 

Li(c^,...,C5^_3)  =  Mi(c^,C5^_g,p\Y') 

(6.29) 

II       II  ti   11      It   II     < 


1  =  ly«.«f 5n-3  • 

By  definition  (4.4) the  complex  quantity 

0"  G   is  a  linear  combination  of  the  coefficients 

n  n  t    t  I  t 

°o'-*"°5n-2*      Pitting     (T^C^  =     <r     +  ir    , 

11    II  II  II 

(T  Cn  =    ^     +  i  T     we  have  relations   of   the  form 

^  -  ^5n-2^'^o'*'*'°5n-2^'  "^  "  ^bn-l^^o' *  *  •'°5n-2^ 
(6.30) 

II  It  "         \         "  _  /    "  "         \ 

^  =  %n-2(®o"**'®5n-2^'    ^  "  ^5n-l^°o' ' '  "^Sn-S^    ' 

At  X  =  -00,  y  =  0   <^(x,y,t)  is  asymptotically 
of  the  form 

$(x,y,t)  =  A  cos(x+t)  +  B  sin(x+t)  +  C  cos(x-t)  +  D  sin(x-t) 
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and  at  X  =  +00  ,  y  =  0  of  the  form 

<^(x,7,t)  =A»cos(x+t)  +  B»sin(x+t)  +  C»cos(x-t)  +  D»sin(x-t)  , 
Here  according  to  (5,18),  {5,19) 

(6.31) 

A'  -  ^-    ■•'  g        *   B'  —  ^'       ' — g 

Prescribing  the  values  of  c  ,Cc__2»^»^»-^*»^* 
we  have  in  (6,29),  (6.30),  (6,31)  lOn+2  eqiiations 
for  the  lOn+2  tuiknomn  quantities 

t      t     "      "     t   t   t   t   11   *i   "   " 

^2.*  *  *  * '^5n-Z>*^±' ** ' *^5n'-Z>*^    '"^   #  O"  »  T  *p  >V  pO"  pZ    • 

If  it  can  be  shown,  that  the  corresponding  homogeneous 
equations  have  only  the  trivial  solution,  then  the 
inhomogeneous  system  will  always  have  a  lonique 
solution.  Thus  there  will  exist  a  unique  motion  of 
the  water  in  the  presence  of  the  barrier  with  energy 
sinks  of  prescribed  strength  at  the  two  points  where 
the  barrier  enters  the  water,  and  corresponding  to 
any  waves  coming  in  from  infinity  on  either  side. 
A  solution  of  the  homogeneous  equations  now 
corresponds  to  a  state  where  ^  ,c  ,Cc_  o^Cgj^  o*-^*» 
B',C,D  vanish.  For  such  a  motion  then  no  logarithmic 
singularities  of  §    exist  {l,e*JJ    (<p^  +  <p^)dxdy  is 
integrable  over  any  finite  region),  and  no  waves  come 
in  from  infinity  from  either  side.  We  have  seen  in 
section  1  that  then  also  A>B,C»,D»  vanish,  and  hence 
<j>'(x,y)  and   ^"(x,y)  vanish  at  infinity.  We  shall 
prove  in  the  next  section  that  then  <Pj^   and   qpg 
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vanish  identically.   This  will  then  complete  the 
proof  of  existence  and  uniqueness  of  a  solution  of 
our  problem  for  given  strength  of  sinks  and  for  given 
incoming  waves.   The  solution  can  be  fo\ind  explicitly 
by  solving  the  lQn+2  linear  equations  indicated  above. 

It  does  not  follow  from  this  argument  that 
Instead  of  the  two  incoming  waves  one  could  prescribe 
the  incoming  and  reflected  wave  on  one  side  at 
infinity,  or  that  there  exists  necessarily  a  standing 
wave  of  given  behavior  at  x  =  -oo  • 
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7.  A  iinlqueness  theorem. 

Theorem:   If  in  addition  to  the  conditions  imposed 

before  the  function  X (z)  is  bounded  everywhere  in 

T  and  vanishes  both  for  x  -*»  +00  and  x  ->  -00  , 

and  if  lim  z  XMz)  =  0>  then  7C  (2)  ^^  constant, 

2  -»•  0 
(It  is  easily  seen  that  the  assumptions  of  this 

theorem  are  satisfied,  if  0   (x,y)  =  R(X)  vanishes 

at  infinity  in  1   ,  and  if   (p(x,y)  is  bounded  near 

0.) 

Proof:   The  assumption  that   X(2)  vanishes  for 
X  =  +00  and  x  =  -00  implies,  according  to  (5,18), 
(5.19)  that 


(7.1)  eC  =  0  ,  Cjj  =  0  . 

Then,  in  conjunction  with  (6,25),  (4.7)  it  is  seen 
that  all  Cj^  vanish  for  k  ■«  0,  ,.,,n.  Lemma  I  on  p.  44 
•  then  guarantees  that  the  functions  B.   {z)   for 
/*•=  0,,,,,n  are  regulaj*  at  infinity.  The  same  holds 
then  for  the  transformed  fvmctions  ^^0  at  the  point     ^ 
Bq.  Using  (7.1),  (6.24),  6.20)  we  have 

(7.2)  Pjj.(w)  =  6Sjj._3^(e^w)  for  k  =  l,...,n 


for  w  in  T-  +  Tn .  Prom  (6,19)  we  have 

(7.3)    P^(w)  =  P   (w)  for   /u  =  0,...,n  . 

As  Pq(w),..,,Pj^(w)  are  regular  at  w  =  B^,  it  follows 
from  (7.2),  that  P^  (w),  .,.,Pj^(w)  are  regular  at  B^, 
and  that  Pq(w),..  .,P  t(w)  are  regular  at  Bg^^^-j^.  Then 
from  (7.3)  Pq(w)  , . ,  .,Pj^(w)  are  regular  both  at  B^  and 
^2n-l'  ®^^  ^^  particular  have  no  branchpoints  there. 


-so- 


Continuing  this  argument  we  find  that  P  {w),,..,P  (w) 
are  regular  at  all  points  B  • 

It  is  known  already  that  the  fxmctions  !E;i,('w) 
are  univalued  and  hence  behave  rationally  near  the 
points  A  •  The  new  assumption  that 

lim  z  %Hz)    =  0 
z->  0 


implies  that 


lim  z^(z)  =  0  . 


Then  G(w)  has  at  most  a  pole  of  order  2n(2n-l)-l  at 
w  =  0  and  a  pole  of  order  at  most  2n-l  at  w  »  oo  •  It 
follows  that  in  the  expansion  (1,29)  of  G(w)  the 
coefficients  c  and  Cc„  o  vanish.  The  vanishing  of 
Ccj.  o  has  as  a  consequence  that  G(w)  is  bounded  at 
infinity.  Then  the  integrands  in  (6,11)  are  of  the 
order  0(-=n),  and  consequently  the  functions  P^(w)  are 

V 

regular  at  w  =  oo  ,  As  those  functions  are  univalued 
at  all  A  ,B  ,00  they  have  no  branchpoint  at  w  =  0 
either.  As  they  then  behave  rationally  at  all 
possible  singular  points,  the  Ptj.(w)  are  rational 
functions.  The  functions  Pij.(w)  were  seen  already  to 
be  rational.  Thus  %  ,  which  is  a  linear  combination 
of  those  functions  is  a  rational  funci^ion  of  w,  as 

X.  is  bounded  in  T"*  it  follows  that  X  is  bounded  in 
T  ,  and  hence  has  no  poles  at  w  =  0,oo,Aq,Bq» 

The  function  K» (w)  has  a  simple  zero  at  each  A  » 
It  follows  that  whenever  X  as  a  fiinction  of  w  is 
regular  at  a  point  A  ,  then  the  two  expressions 

(7.4)     eX»  =  £— ^Is  a^d  f  =  X^  +  iA 


-81- 


liave  a.t  most  poles  of  first  order  at  that  A  . 
Inversely,  v^henever  one  of  the  two  expressions 
(7,4)  has  a  pole  of  first  order  at  a  point  A^, 
then  X  is  regular  at  that  point.  Now  6X ' 
is  real  along  the  line  0,Aq,oo,  and  hence  has  the 
same  behavior  at  the  points  A  and  Ao  ^^,   It  follows 
that  if  X  is  regiilar  at  A^,  then  also  at  Agj^^p* 
Similarly  the  function  f  =  K,^    +  ±X     Is  real  along 
the  positive  real  w-axis,  and  hence  has  the  same 
behavior  at  the  points  A  and  Ag^^^^^-j^,  Hence,  if  X, 
is  regular  at  A  ,  it  is  regular  at  Ag^^^^^-j^.  Hence, 
as  X  is  regular  at  A  ,  it  is  regular  at  all  A^. 
It  follows  that  the  rational  function  X  ^^s  no 
poles  at  all  in  the  w-plane,  and  hence  is  a  con- 
stant. This  completes  the  proof  of  the  theorem. 
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